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ABSTRACT 

In this paper we use the 7 mm and 1.3 mm light curves obtained during the 
2003.5 low excitation phase of the rj Carinae system to constrain the possible 
parameters of the binary orbit. To do that we assumed that the mm wave 
emission is produced in a dense disk surrounding the binary system; during 
the low excitation phase, which occurs close to periastron, the number of 
ionizing photons decreases, producing the dip in the radio emission. On the 
other hand, due to the large eccentricity, the density of the shock region at 
periastron is very high and the plasma is optically thick for free-free radiation 
at 7 mm, explaining the sharp peak that was observed at this frequency and 
lasted for about 10 days. From the shape and duration of the peak we were able 
to determine the orbital parameters of the binary system, independently of 
the stellar parameters, such as mass loss rates, wind velocities or temperature 
at the post-shock region. 

Key words: stars: individual {rj Car) binaries: general stars: variable radio 
continuum: general 



1 INTRODUCTION 

The occurrence of the predicted 2003.5 low excitation event in the rj Carinae system, ex- 
tending from radio to X-rays left little doubt about its periodic behavior (Lajiis et al. 2003, 
van Genderen & Sterken 2003, Smith et al. 2004, Whitelock et al. 2004, Abraham et al. 
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2005, Corcoran 2005, Stahl et al. 2005, Weis et al. 2005). Although the binary nature of 
the system seems to explain the strict periodicity as well as the high intensity of the X-ray 
emission, attributed to wind- wind collision (Pittard et al. 1998, Ishibashi et al. 1999), orbital 
parameters, like eccentricity, epoch of periastron passage and orientation of the orbit rela- 
tive to the observer are not well defined (Davidson 1997, van Genderen et al. 1999, Pittard 
& Corcoran 2002). Besides, other observed properties like the optical and UV light curves 
and the recently discovered He II A 4686 A emission (Martin et al. 2004, Steiner & Damineli 
2004) seem to be more appropriately described by periodic shell events. 

Falceta-Gongalves, Jatenco-Pereira & Abraham (2005, Paper I) were able to conciliate 
the two hypothesis, still assuming a binary system to explain the strong X-ray emission, but 
also taking into account that, near periastron and because of the highly eccentric orbit, the 
wind behind the shock cools rapidly and forms dust grains, which absorb part of the optical 
and UV flux, as in a shell-ejection event. For this process to be effective, at periastron the 
secondary star must be located between 77 Carinae and the observer, solving also the dis- 
crepancy between the orbital parameters derived from ground and space based observations 
of the radial velocities of the H recombination lines (Davidson 1997). 

Although radio waves are not affected by dust absorption, the radio light curves show 
also periodic dips, similar to those found at X-rays (Cox et al. 1995, Abraham & Damineli 
1999, Duncan, White & Lim 1997, Abraham et al. 1999). Based on high resolution radio 
observations, Duncan & White (2003) suggested that the radio emission is produced by the 
free-free process in an extended (2" radius) disk surrounding the r\ Carinae system; the dips 
would be the consequence of a sudden decrease in the available ionizing UV flux, as expected 
in a shell-like event. 

To better understand the mm-wave light curve during the 2003.5 event, Abraham et 
al. (2005) made daily observations at 7 mm with the Brazilian Itapetinga radiotelescope 
and weekly observations at 1.3 mm with the ESO (European Southern Observatory) SEST 
radiotelescope at La Silla, Chile. The closely spaced observations provided not only the 
light curves with great precision, but also showed an unexpected increase in flux density, 
which could not be explained by free-free wind emission from the individual stars. Using 
reasonable values for the physical conditions at the shock during periastron passage, they 
showed that the observed excess flux density is compatible with optically thick free-free 
emission produced at the wind collision site. 

In this paper we show that the mm-wave light curve during periastron passage, after 
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subtracting the disk contribution, can be explained by simple geometric arguments, if the 
free-free emission originates in the conic surface defined by the wind-wind collision interface. 
Furthermore, we show that the orbital parameters are strongly bound by the light curve, 
and that periastron should occur when the secondary star is between r\ Carinae and the 
observer, as also required to explain the X-ray light curve and shell-like events (Paper I). 

In Section 2 we describe the emission model and the geometry of the shock surface. In 
Section 3 we present the results including expressions for the emission at different positions 
in the secondary's orbit and the determination of the orbital parameters. In Section 4 we 
discuss our results and compare them with those commonly accepted for the r\ Carinae 
system and in Section 5 we summarize our conclusion. 



2 MODEL FOR THE MM- WAVE EMISSION 

The observed light curves, obtained with the Itapetinga radiotelescope at 7 mm and with 
SEST at 1.3 mm are presented in Figure 1, where we can see the decrease in flux density 
superposed to a sharp peak at 7 mm, which reaches its maximum intensity on 29 June 2003. 
Details of the observations are found in Abraham et al. (2005). 

We will assume that the observed free-free emission is formed by the contribution of 
both an extended ionized disk, which surrounds the binary system, and the wind-wind 
shock material. At the beginning of the low excitation phase there is a sudden decrease in 
the UV ionizing flux, as proposed by Duncan & White (2003); the central denser part of the 
disk recombines in timescales of a few days and the flux density decreases. Also, since the 
recombination timescale for the disk material depends on the unknown density distribution 
and we are interested mainly in the shock emission, for which the physical conditions are 
constrained by the orbital parameters, we will model the disk emission by an exponential 
decay function of the form: 

Sdisk(7mm) = a exp(bt + c) + mt + n (1) 

Abraham et al. (2005) determined the free parameters a, b, c, m and n by fitting the 
data at both sides of the peak in the 7 mm light curve. They also calculated the physical 
conditions at the shock site required to obtain the observed peak flux density of about 1 Jy 
at this wavelength. The most favorable situation involved a high temperature optically thick 
plasma, as expected to be produced during periastron passage. In this paper we will obtain 
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Figure 1. mm- wave light curves of r\ Carinae at the epoch of the predicted low excitation phase. Open circles correspond the 
1.3 mm emission (right axis), full circles correspond to 7 mm emission (left axis). 
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Figure 2. Surface of momentum balance between too colliding winds in a binary system. The stars are indicated by dots. The 
separation between them corresponds to the 77 Carinae system at periastron passage for an eccentricity e = 0.95. The aperture 
angles of the asymptotic cone are 56° and 44° for r) values of 0.2 and 0.1, respectively. 



simultaneously the orbital and disk emission parameters that reproduce the total observed 
light curves, both at 1.3 and 7 mm. 



2.1 The free-free shock emission 

At mm wavelengths, the flux density due to free-free emission must be calculated taking 
into account radiation transfer since, at the densities and temperatures found in the shock, 
the optical depth can be large. Using the Rayleigh- Jeans approximation we have: 

dS (\) = ^(l-e~^)dn, (2) 
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Figure 3. Left: schematic view of the binary orbit near periastron. The shock region is shown as a projected cone: the 
secondary and r] Carinae shocks are represented by the dark and the hatched regions, respectively. Points (2) and (4) represent 
the positions in the orbit for which the line joining the stars is perpendicular to the line of sight. The position of the observer 
at conjunction or opposition is represented by point (3). The angle (f> represents the position angle of the secondary in its orbit, 
with r] Carinae at the focus of the ellipse. Periastron passage corresponds to <f> = 0°. Right: artist view of the emitting surface 
with $ = 0. 

where dfl is the element of solid angle subtended by the emitting region and r(A) the optical 
depth, calculated from: 

r(A)= [ L K ff (X)ds, (3) 

where L is the depth of the emitting region and ft//(A) is the free-free absorption coefficient, 

which at radio wavelengths can be calculated from: 

3.7 x 10 8 h n t n e ^ // (A)A 2 
K /A A ) = ~ 2 ^372" ' W 

and <?//(A) is the Gaunt factor for free- free emission. For the physical conditions at the 
shock, T > 10 5 Z 2 and hc/X <C kT, 57/ (A) can be approximated by: 

The optical depths at 1.3 and 7 mm can be determined from the ratio of their respective 
flux densities. Although the sharp peak in the light curve was not observed at 1.3 mm, 
from Figure 1 we can estimate the contribution of the shock emission at this wavelength. 
The emission from the extended disk is about 5 Jy at the minimum and the total 1.3 mm 
flux densities before and after the 7 mm maximum are about 11 Jy and 9 Jy, respectively. 
The difference gives the approximate value of 5 Jy for the expected shock contribution. 
Assuming a flux density of 1 Jy for the 7 mm peak emission, we find r(7 mm) ~ 4 and 
r(1.3 mm) ~ 0.2. Therefore, close to the peak, the 7 mm emission is optically thick and the 
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form of the light curve is completely determined by the variation of solid angle subtended 
by the shock surface at each position of the binary orbit. 

2.2 The geometry of the shock surface 

The general form of the momentum balance surface between the two colliding winds was 
given by Stevens, Blondin & Pollock (1992) as: 

dy (y- 1/2 d 2 2 + d x 2 )y 

dz r 1 - 1 lH 2 2 z + d 1 2 (z- D)' 1 ' 

where D is the distance between the stars; d\ and d 2 are the distances of the primary and 
secondary star to the contact surface, respectively; the coordinates z and y are measured 
along the line that joins the stars and perpendicular to it, respectively; rj = M s v s /M p v p , 
where M p and M s are the mass loss rates of rj Carinae and the companion star, v p and v s 
their respective wind velocities. 

The surface can be described asymptotically by a cone with an opening angle j3 that 
depends on rj. Figure 2 shows a schematic view of the contact surface for two values of rj 
and a separation between the two stars equivalent of what is expected at periastron in the r\ 
Carinae binary system. The asymptotic angle defined by the surface is 56° for rj = 0.2 and 
44° for rj = 0.1. 

When radiative losses are high, a thin and dense shock layer will form at both sides and 
very close to the contact surface (Eichler & Usov 1993). Due to the orbital motion, the cone 
will be distorted beyond a distance from the cone origin r > (P/2)v p , where P is the orbital 
period. For the rj Carinae system, with a period of 5.5 years and wind velocity of 700 km s _1 , 
rsi6x 10 15 cm. 

In Paper I we assumed M p = 2.5 x 10~ 3 M Q yr _1 , v p = 700 km s _1 for rj Carinae and 
M s = 10~ 4 M yr _1 , v s = 3000 km s -1 for the secondary star, which results in rj = 0.2. 
We calculated the temperature and density at both the primary and secondary shocks near 
periastron as T p = 6 x 10 6 K, n p = 2 x 10 11 cm~ 3 and T s = 2 x 10 8 K, n p = 5.6 x 10 9 
cm -3 , respectively. Also, we pointed out that the actual temperature can be lower and, if we 
assume isobaric equilibrium, the density higher. Using these values, Abraham et al. (2005) 
estimated a radius for the projected surface necessary to reproduce the 7 mm peak emission 
of about 10 14 — 10 15 cm. Under these circumstances, the shock surface can be approximated 
by the non distorted asymptotic cone. 
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In Figure 3 we present a schematic view of the contact surface position and orientation 
with respect to the observer at different orbital phases, scaled for eccentricity e = 0.9. is the 
position angle of the secondary star orbit relative to rj Carinae, with = at periastron and 
= $ at conjunction. The dark and hatched regions represent the secondary and primary 
shocks, respectively. As we can see, if periastron occurs close to conjunction, for | 0+$ | < (3 
only the secondary shock is visible to the observer, for 90° > | + $ | > (3 only part of 
the emission from the secondary shock is directly visible, the remaining part is partially 
absorbed by the primary shock, which also contributes with its own emission. For all other 
position angles, only the primary shock surface is directly observed. If periastron occurs close 
to opposition, the situation is reversed with respect to the secondary and primary emission. 

3 RESULTS 

3.1 Shock emission along the orbital phase 

For orbital inclination 90° ^ % > /3, the observed emission is produced in the internal and 
external parts of the contact cone, which is located at a distance D from the observer, and 
has an extension R{4>) ■ The optical depth of the primary and secondary shock regions are not 
known, except for the fact the maximum in the 7 mm light curve requires r s ~ 4 (Abraham 
et al. 2005). We will assume r SjP oc [c? s , p (0)/(i s ,p(O)]" a! , where rf SjP (0) represents the distance 
from the shock to the respective star, with r s (0) = 4; a and r p (0) are free parameters to be 
determined from the fitting to the observed light curve. We will also assume that the optical 
depths are constant along the solid angle subtended by the emitting surface (in fact t SjP (0) 
is the average over the solid angle). We will also assume that R(4>) varies as [c4 iP (0)/c4 iP (O)] 2 
as estimated by Eichler & Usov (1993) and % = 90°, as seems to be implied by the high 
resolution radio observations of Duncan, White & Lim (1997). Smaller inclination angles, 
in the range 90° ^ i > {3 will affect the projected surface area, requiring a larger value of 
R, which is anyway a free parameter in our model. We present the expressions for the total 
emission as a function of the orbital phase (— n ^ ^ it) for the case in which periastron 
is close to conjunction, for the opposite case the role of the secondary and primary shock 
emission should be exchanged. We include expressions for the emission from the primary 
shock, although its contribution is negligible compared to that of the secondary shock when 
the latter arises from an optically thick region, since T s 3> T p . 

For | + $ | < 13: 
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^(0) = ^[l- e - rsW ]^i(0), (7) 
with 

= irab(<f)) (8) 
where a = (R/D) sin (3 and 5(0) = (R/D) sin/3 | cos(0 + $) |. 
For 90° > | + $ | > (3 

S(0) = S ls (0)+S lp (0), (9) 
with 

IkT 

SM = -^[1 - e-^)][Q 2 (0) + e^ W fi 3 (0)], (10) 
and 

9 FT 

5i P (0) = ^[l-e-^W]Q 3 (0), (11) 
where 

ft 2 (0) = a6(0)[ tt/2 + (1 - C 2 ) 1/2 + arcsinC ], (12) 
and 

fi 3 (0) = a fe(0) S in /3 cos (3 (1 - C 2 ) 1/2 [ | sin(0 + $) | -( 2 ], (13) 
with ( = | tan/3/tan(0 + <E>) |. 
For 180° - (3 ^ | + $ | > 90°: 

S(0) = S 2s (0) + S 2p (0), (14) 
where 

S 2s (0) = -j^e-^[l - e-^mM + Q 3 (0)], (15) 
and 

5 2p (0) = -^[1 - e^ W ][fi 2 (0) + n 3 (^)], (16) 
For | + $ |> 180° 

S(0) = S 3s (0) + S 3p (0), (17) 
where 

^3s(0) = ^e-^[l - e^W]Q(0), (18) 



Radio emission from the colliding winds of r\ Carinae 9 



3.0 



2.5 - 



„ 2.0 



E 1.5 
E 

w 1.0 




0.5 - 




20 

time (days) 




-2 1 — ; — ' — — ^ — 1 1 1 1 ' 1 -1 

-20 20 40 60 

time (days) 



Figure 4. Model of the 7 mm light curve during the low excitation phase. Upper panel: fitted model together with the observed 
7 mm flux density (circles); the disk model emission is also shown. Lower panel: shock emission model (continuous line) and the 
observed shock contribution (circles), calculated as the difference between the total observed flux density and the disk model 
contribution. The residuals of the fitting, calculated as the difference between the observation and the total model emission 
(disk plus shock) are also shown. Vertical lines represents orbital position angles from left to right: + = —90°, — /?, j3 and 
90°. Zero in the time axis corresponds to June 29, 2003 (JD = 2452819). 



and 

S 3p (0) = ^[l-e-^M0) 



(19) 



3.2 Determination of the orbital parameters from the mm-wave light curves 

We fitted the observed 7 mm light curve simultaneously to the shock emission model de- 
scribed in the previous subsection, using the parameters listed in the first column of Table 
1, and to the surrounding disk emission, described at the beginning of Section 2. We found 
that each parameter affects differently the form of the light curve, in such a way that they 
do not let too much margin for variation. 

The width and shape of the light curve near the maximum were affected mainly by the 
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orbital parameters: the eccentricity e determines how fast the secondary's position angle 
changes with time, the phase of conjunction $ determines the asymmetry in the shape of 
the peak and the epoch of periastron passage defines its position in the light curve. 

The aperture angle of the emitting cone (3 together with the rate at which the opacity 
changes, given by the parameter a, were responsible for the shape of the light curve at larger 
phase angles. For orbital angles | + $ | > (3, the form of the light curve is also affected 
by the absorption and emission of the primary shock; if these last effects were important we 
should see a discontinuity in the light curve at the epoch in which | <fi + $ | = /3, unless the 
absorption of the secondary shock emission by the primary is exactly compensated by the 
emission from the primary shock. Since we did not see any discontinuity and we know that 
the emission from the primary shock is much smaller than that of the secondary, because 
its temperature is much lower, we neglected both absorption and emission from the primary 
shock. 

The parameter range for the model shock emission that allowed us to reproduce the 
observed light curve are presented in the second column of Table 1, the parameters of the 
best model are listed in the third column. The best parameters for the disk emission model 
represented by equation (1) are: a = 3.1 x 10~ 5 , b = -0.228, c = 7.50, m = -0.004 and 
n = 0.94, the time t was measured from 29 June 2003 (JD = 2452819). 

We found the best fitting for e = 0.95 and a = 4, although solutions were also possible 
for 0.93 ^ e ^ 0.95, in which case 5.5 ^ a ^ 4. For the best model the orbital angle 
at periastron is $ = —30° ± 5°, while for the other eccentricities it varied in the interval 
—50° ^ $ ^ —30°. For such high eccentricities the differences in the epoch of periastron 
passage and conjunction relative to the maximum in the light curve (June 29, defined as 
t = in our graphs) are very small: 3.2 ^ t(4> = 0) ^ 1.5 and —1 ^ t(<j) = $) ^ 0.1 
respectively, where t is measured in days. 

The comparison of the model with the observations can be seen in Figures 4. The upper 
panel shows the total fitted model emission (disk plus shock), together with the observed 
7 mm flux density (circles) and the model disk emission. The lower panel shows the shock 
emission model (continuous line) and the observed shock contribution, calculated as the 
difference between the total observed flux density and the disk model contribution. This 
panel also shows the residuals, calculated as the difference between the observed flux density 
and the added disk and shock model contribution. 

Although the 1.3 mm light curve has not enough temporal resolution to show a peak in 
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Figure 5. Model for the 1.3 mm light curve. Upper panel: 1.3 mm data (filled circles) and model shock emission (solid line), 
obtained with the parameters presented in Table 1; lower panel: expected disk emission, both at 1.3 and 7 mm, calculated as 
the difference between the observations and the shock model flux density. Zero in the time axis corresponds to June 29, 2003 
(JD = 2452819). 



the light curve, we know that the predicted emission must be proportional to the optical 
depth at each point in the orbit since, as discussed before, the region should be optically 
thin at this wavelengths. As in the 7 mm model, we assumed that the optical depth varies as 
[d StP (4>) I d StP (0)]~ a , its value at periastron can be adjusted to give the right flux density and 
a is already determined by the 7 mm emission model. The model emission at 1.3 mm for 
a = 4 can be seen in the upper panel of Figure 5, together with the observations. The lower 
panel shows the expected disk emission, both at 1.3 and 7 mm, calculated as the difference 
between the observations and the shock model flux density. 
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Table 1. Parameters for the r) Carinae binary system derived from mm- wave shock emission 
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Figure 6. Model flux densities for eccentricities 0.9 (light line) and 0.95 (heavy line), together with the observed flux density 
(circles). Zero in the time axis corresponds to June 29, 2003 (JD = 2452819). 

4 DISCUSSION 

The orbital parameters of the t] Carinae binary system were estimated under the assumption 
that the mm-wave emission during the low excitation phase is produced both by an ionized 
disk surrounding the system and by the shock region formed by wind-wind collision. 

In our model for the shock emission the eccentricity is the best defined parameter; small 
differences in its value cause a large broadening of the radio light curve. As mentioned in 
the last section, fitting to the light curve were obtained for 0.93 ^ e ^ 0.95, these values are 
slightly higher than the value of 0.9 found by Pittard & Corcoran (2002) in their numerical 
simulations of the X-ray spectral behavior; however this small difference produces large 
changes in our 7 mm shock emission model, as can be seen in Figure 6. These differences 
cannot be compensated by a decrease in the cone aperture angle f3 or an increase in the 
dependence of the optical depth with the distance between the shock and the star. The 
aperture angle of the cone f3, which represents the shock surface asymptotic behavior, was 
varied between 40° and 60°. These values correspond to 0.1 < rj < 0.2, similar to those 
found by Pittard & Corcoran (2002). Due to the high eccentricity, the epochs of periastron 
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passage and opposition are very close, differing in only a few days, even though the phase 
for opposition is not small (—30° < $ < —50°). 

The fitting of the light curve described in the last section was made under the assumption 
that periastron occurs close to conjunction, that is, when the secondary star passes between 
the observer and rj Carinae. If periastron occurs close to opposition, the role of primary and 
secondary shock emission are reversed. The behavior of the light curve for | 0+$ | ^ (3 would 
be similar, but for larger phase angles the contribution of the secondary shock emission, 
much higher than that of the primary, would produce an increase in the model light curve, 
incompatible with the observations. We consider this result as additional evidence for the 
case presented in Paper I, in which we explained the X-ray light curve and shell-like behavior 
of the r] Carinae binary system by assuming the secondary star positioned between r] Carinae 
and the observer close to periastron passage. 

The eccentricity values found in our model are compatible with the variation of the 
H recombination lines radial velocities with orbital phase, as measured by Damineli et al. 
(1997, 2000) and Davidson (1997), if we assume that the lines are produced in the primary 
shock instead of in the stellar atmosphere (Hill, Moffat & St-Louis 2000, 2002). With these 
large values for the eccentricity, the distance from the shock to rj Carinae at periastron turns 
out to be less than 1 AU and, depending on the extension of its atmosphere, it can change 
the shock behavior on the primary side. However, we do not expect too much changes in 
our emission model, which depends mainly on the secondary shock. 

Finally, the linear extension R of the shock cone can be calculate from the actual value 
of the observed 7 mm flux density, resulting in R = 4 x 10 14 cm for a plasma temperature 
of 10 7 K. This value is valid at the orbital phase corresponding to the peak of the 7 mm 
emission. If this value were not constant along the orbit, the width of the peak would be 
different. In the analytical treatment of the shock properties, Usov (1991) found that the 
size increases with the distance between the stars, which will broadened the peak model. 
However, for | + $ | < (3 the difference in distance is less than a factor of two 

5 CONCLUSIONS 

In this paper we used the 1.3 and 7 mm light curves of rj Carinae during the 2003.5 low 
excitation phase to determine the orbital parameters of the binary system. We assumed that 
the emission is produced both in the extended disk surrounding the 77 Carinae binary system, 
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as proposed by Duncan & White (2003), and at the shock region, the later responsible for 
the sharp peak seen in the 7 mm light curve during the decreasing phase of the disk emission. 
We calculated the form of the momentum equilibrium surface and the physical conditions 
in the material behind the primary and secondary shocks and showed that the emission 
could be produced in the non distorted part of the asymptotic cone formed by the two 
shocks, as proposed by Abraham et al. (2005). From the ratio of the 7 mm and 1.3 mm 
observed flux densities, we concluded that the emitting region should be optically thick at 
7 mm and optically thin at 1.3 mm. For this reason, the shape of the 7 mm light curve 
close to periastron was completely determined by the value of the solid angle subtended by 
the emitting surface, which changes with the orbital position angle of the secondary star, 
and therefore, depends only on the orbital parameters. We found satisfactory agreement 
between models and observations for orbits with eccentricities 0.93 ^ e ^ 0.95, orbital angle 
at periastron —50° ^ $ ^ —30° and cone aperture angle 40° < (3 < 60°, compatible with 
what was expected from other observations, like X-ray fluxes and radial velocities of the H 
recombination lines (Pittard & Corcoran 2002, Damineli et al. 2000, Davidson et al. 2000). 

The epoch of periastron passage and opposition varied between 3.2 ^ t(<fi — 0) ^ 1.5 
and —1 < t(4> = $) ^0.1 respectively, where t is measured in days from June 29, 2003. 

Finally we must mention that the model presented here is a very simplified version of a 
very complex process and that confirmation of the orbital parameters presented here should 
be obtained from detailed numerical simulations. 

ACKNOWLEDGMENTS 

This work was supported by the Brazilian Agencies FAPESP, CNPq and FINEP. 
REFERENCES 

Abraham Z., Damineli A., 1999, in Carinae at the Millenium, ASP Conference Series 179, 
San Francisco, p263 

Abraham Z., Damineli A., Durouchoux, P., Nyman L-A., McAuliffe F., 2002, in Cosmic 
Masers: from Proto-stars to Black Holes, ASP Conference Series 206, San Francisco, p234 
Abraham, Z., Falceta-Gongalves, D., Dominici, T. P. et al. 2005, A&A, 437, 977 
Churchwell, E., Bieging, J. H., van der Hucht, K. A. et al., 1992, ApJ, 393, 329 
Corcoran, M. F. 2005, AJ, 129, 2018 



Radio emission from the colliding winds of i] Carinae 15 

Cox P., Mezger P. G., Sievers A. et al., 1995, A&A, 297, 168 179, 
Damineli, A. & Conti, P., Lopes, D., 1997, NewA, 2, 107 

Damineli, A., Kaufer, A., Wolf, B., Stahl, O., Lopes, D. & Araujo, F. 2000, ApJ, 528, L101 
Davidson K., 1997, NewA, 2, 397 Astrophys., 35, 1 

Davidson, K., Ishibashi, K., Gull, T. R., Humphreys, R., Smith, N., 2000, ApJ, 530, L107 
A., 319, 2003, 

Duncan, R. A., White, S. M., Lim, J., 1997, MNRAS, 290, 680 
Duncan, R. A., White, S. M. 2003, MNRAS, 338, 425 
Eichler, D., Usov, V. 1993, ApJ 402, 271 

Falceta-Gongalves D.,Jatenco-Pereira V., Abraham Z, 2005, (Paper I) MNRAS, 357, 895 

Hill G., Moffat A., St-Louis A., 2000, MNRAS, 318, 402 

Hill G., Moffat A., St-Louis A., 2002, MNRAS, 335, 1069 ApJ, 

Ishibashi K. Corcoran M., Davidson K., Swank J., Petre R., Drake S., Damineli A., White 

S., 1999, ApJ, 524, 983 
Lajiis, F. E., Gamen, R., Schwartz, M. et al. 2003, IBVS, 5477, 1 
Martin, J. C, Koppelman, M. D. et al. 2004, ApJ, 127, 2352 D. 
Pittard J. M., Stevens I. R., Corcoran M. F., Ishibashi K., 1998, MNRAS, 299, L5 
Pittard J. M., Corcoran M. F., 2002, A&A, 383, 636 

Smith, N., Morse, J. A., Collins, N. R., Gull, T. R. 2004, ApJ, 610, L105 medium, 

Stahl O., Weis, K, Bomans, D. J. et al. 2005, A&A, 435, 303 

Steiner, J. E., Damineli, A. 2004, ApJ, 612, L133 

Stevens, I. R., Blondin, J. M., Pollock, A. M. T., 1992, ApJ, 386, 265 

Usov V. V., 1991, MNRAS, 252, 49 

van Genderen, A. M., Sterken, C, de Groot, M., Burki, G. 1999, A&A, 343, 847 
van Genderen, A. M., Sterken, C. 2003, A&A, 423, LI 
Weis, K., Stahl, O., Bomans, D. J. et al. 2005, AJ, 129, 1694 
Whitelock, P. A., Marang, F., Crause, L, Corcoran, M., 2003 IAUC, 8160 
Whitelock, P. A.,Feast, M. F., Marang, F., Breedt, E. 2004, MNRAS, 352, 447 



